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Introduction 

Consider a weight on R, depending on t = (ti, £2, • • •) G C 
p t {z)dz = e£ ttzi p(z)dz = e - y(z)+ Sr Wdz, with 



m v , (z) 9(Z 



/(*)' 

(0.1) 



Hermitean matrix integrals (revisited) This weight leads to a t-dependent 
moment matrix 



m n (t) = {Hk+t{t))o<k,e<n-i = (^J^z k+e p t (z)dz^ 



0<kl<n-l 



with the semi- infinite moment matrix m^, satisfying the commuting equa- 
tions 

dm^/dtt = AVoo = m^A*, (0.2) 

where A is the customary shift matrix. Considering the lower- and an upper- 
triangular matrix Borel decomposition 



moo = S S 



1 oT-1 



(0.3) 



which is determined by the following t-dependent matrix integrals^] (n > 0) 

T„(t) := / e Tr( - y(x)+ ^' lX!) rfX = detm n , and r = 1, (0.4) 

with Haar measure dX on the ensemble 7i n = {n x n Hermitean matrices}. 
As is well known, the integral (0.4) is a solution to the following two systems, 

(i) the KP-hierarchy 

( ~ I d 2 \ 

ip k+4 (d) - 2 dtidt ) T n°r n = 0, for k, n = 0,1,2,.... (0.5) 

(ii) the Toda lattice; i.e., the tridiagonal matrix 

1/2 



(A 

dt 



L(t) := SAS- 1 



- log ^ (^) 

1 ° TO V T( J 

1/2 » 

ati <=> n 



~l ) ~1 



1/2 







N 1/2 

T1T3 \ 



OTl T2 



V 



(0.6) 



'We set vol(W(n)) = 1 for all n. 
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satisfies the following commuting Toda equations 



dL_ 



where (A) s k denotes the skew-part of the matrix A for the Lie algebra split- 
ting into skew and lower-triangular matrices. Moreover, the following t- 
dependent polynomials in z, are defined by the S-matrix obtained from the 
Borel decomposition (0.3); it is also given, on the one hand, by a classic 
determinantal formula, and on the other hand, in terms of the functions 
T n (t): 



( 



p n (t,z) : = (S(t)x(z)) r 



\J T n T n+l 



det 



m n (t) 



(J>n,n-l{t) 



1 \ 

Z 



z n ) 



,Jn{t-[z- 1 ]) 



\/ T n T n+l 

The p n J s are orthonormal with respect to the (symmetric) inner-product 
(,) sy , defined by (z l ,z^) sy = pij, which is a restatement of the Borel de- 
composition (0.3). The vector p(t, z) = (p n (t,z))n>o is an eigenvector of the 
matrix L(t) in (0.6): 

L(t)p(t,z) = zp(t,z). 

Symmetric and symplectic matrix integrals. Instead consider the 
following skew-symmetric matrix = (^ij)ij>o of moments^ 



a4?(*) 



(2) 

x l y 3 e(x-y)p t {x)p t {y)dxdy or p>J{t) 



R 



{y\y j }pt(y) 2 dy, 

(0.7) 



both satisfying the equations 

drrioo 



A'rrioo + m^A 



Ti 



In this paper, we consider symmetric matrix integrals^ 

r£(t) := J- [ P ("VW+Er^dX = pf(m£), 
(2n)! Js 2n 



(0.8) 



(0.9) 



2 e(x) = 1, for x > and = — 1, for x < and {/, g} = f'g - fg'. 

3 where again we set the volume of the orthogonal and symplectic groups equal to 1. 
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and symplectic matrix integrals 

T&\t) := i / e 2 ^ l-vW+XWdX = p/(mg), (0.10) 

^' •'Tin 

both expressed in terms of the Pfaffian of the "moment" matrix m^, where 

(1) in the first case, dX denotes Haar measure on the space S 2n of symmetric 
matrices and, 

(2) in the second case, dX denotes Haar measure on the 2n x 2n matrix real- 
ization T 2n of the space of self-dual n x n Hermitean matrices, with quater- 
nionic entries. 

Since moo is skew- symmetric, the Borel decomposition of moo will require the 
interjection of a skew-symmetric matrix J, used throughout this paper, 



J 



\ 



1 
-1 



1 
-1 







1 
-1 



V 



with J = -I 



(0.11) 



/ 



and the order 2 involution on the space V of infinite matrices 

J : V — ► V : a i — > J(a) := J a 1 J. 
The skew-Borel decomposition 

moo(t) = Q~ l (t)JQ~ 1T (t), 



(0.12) 



(0.13) 



can entirely be expressed in terms of the integrals T2 n (t), (0.9) and (0.10) 
corresponding respectively to the first and second moment matrix (0.7). They 
satisfy both 

(i) the Pfaffian KP-hierarchy for k, n — 0, 1, 2, . . . , 



1 d 2 ' 
2ati9t fc+ 3 



(0.14) 
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(ii) the Pfaff lattice; i.e., the matrix, constructed by dressing up A with Q 
and which this time is full below the main diagonal, 



L = QAQ- 1 = h- 1/2 



( ^ 00 














Ll2 





* 


L>21 


L22 


L23 


* 


* 


L32 


L33 



V 



til\ 



satisfies the Hamiltonian commuting equations 



dL 

dti 



({D) + + J{{L% 



+ l((L% + J((L%)),L 



with the entries and the entries of h, being 2x2 matrices 
h = di&g(h I 2 , h 2 h, hj 2 , ■ ■ ■), h 2n = T 2n+2 /T 2n , 



-(logr 2n )- 



S2(d)T 2n _ S'2(-9)T2n+2 
T2n T2n + 2 



1 \ 

(l0gT 2 „ +2 )7 



J n,n+1 



(0.15) 



0" 

1 0. 



J n+l,n • 



* (logr 2n+2 )- 



(0.16) 



The following t-dependent polynomials q n (t,z) = (Sx(z)) n in z, defined by 
the S matrix of the skew-Borel decomposition (0.13), have determinantal and 
Pfafnan r-function expression, in analogy with the Hermitean case: 



<?2n(M) 



( 



1 



pf 



= z 



V / r 2n {t)T2n+2(t) 
2n T 2 „(t- [Z- 1 ]) 



-1 



m 2 n+i(t) 



—Z 



2ri 



1 \ 

Z 

^2n 
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and 



Q2n+i(t, z) 



V /r 2n r 2n+2 



Pf 





1 


^0,271+1 




Z 


/^l,2n+l 


m 2n (t) 


z 2n 


/^2n,2n+l 


-1 -z 





_~2n+l 


~/^0,2n+l —y-X,2n+l ■ ■ ■ 


z 2n+l 






<9 



^r 2n {t)r 2n+2 {t) \ dtij 
y 1 r 2n {t)T 2n+2 {t) 



(0.17) 



form skew-orthonormal sequences with respect to the skew inner-product 
(, ) s k, defined by (y l ,z^) s k = fiij, namely, we have the following restatement 
of the skew-Borel dcomposition (0.13): 

{(qi,Qj))o<i,j<oo = J- (0.18) 

Finally, the vector q(z) = (q n (z)) n >o forms a eigenvector for the matrix L: 

L(t)q(t,z) = zq(t,z). (0.19) 

In section 2, we show how a general skew-symmetric infinite matrix flow- 
ing via (0.8) and its skew-Borel decomposition (0.13), lead to wave vectors 

satisfying bilinear relations and differential equations. Section 3 deals 
with the existence, in the above general setting, of a so-called Pfaffian r- 
function, satisfying bilinear equations and a KP- type hierarchy. In |J, these 
results were obtained, by embedding the system in 2-Toda theory, while in 
this paper, they are obtained in an intrinsic fashion. 

For k = 0, the KP-like equation (0.14) has already appeared in the context 



of the charged BKP hierarchy, studied by V. Kac and van de Leur |fL2| ; the 
precise relationship between the charged BKP hierarchy of Kac and van de 
Leur and the Pfaff Lattice, introduced here, deserves further investigation. 
See the recent paper of van de Leur |L6|. 



A remarkable map from Toda to Pfaff lattice: Remembering the no- 
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tation (0.1), we act with the ^-operator, 

on the t-dependent orthonormal polynomials p n {t, z) in z; in @, we showed 
that the matrix M defined by 

n t p(t, z) = (f\L)M - ^-(L)^ p(t, z) =: Np(t, z) (0.21) 

is skew-symmetric. The t-dependent matrix M is expressed in terms of L 
and a new matrix M, defined by 

zp = Lp and e ~^ tkzk -^-e^ tkzk p = Mp. (0.22) 

dz 

Consider now the skew-Borel decomposition of Af(2t) and its inverse^] Af(2t)~ 1 , 
in terms of lower-triangular matrices Ot + \ (t) and Ot-\ (t) respectively^] : 

A/^ 1 = -O^JOj^it). (0.23) 

Then, the lower-triangular matrices 0(±)(i) map orthonormal into skw- 
orthonormal polynomials, and the tridiagonal L-matrix into an L-matrix: 

p n (t,z) .— > q^\t,z)=O i±) (t)p n (t,z) 

L(t) .— > L(t) = O i±) (t)L(2t)O w (t)- 1 . (0.24) 

(Toda Lattice) (Pfaff Lattice) 

It also maps the weight into a new weight 

= e -vw ^ = e -vW := e -i(vW T iog/W) j 

and the corresponding string of r-functions into a new string of pfaffian r- 
functions: (remember V t (z) = V(z) — U^) 



rt\t):= / e^W>to, 09 = 4) 

^ r«(t):= / e T ^~ v ^dX ((3=1). 

4 See the appendix, section 9. 

5 The upper-signs (respectively, lower-signs) correspond to each other throughout this 
section. 
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For the classical orthogonal polynomials p n (z), we have shown in 
that jV(0) is not only skew-symmetric, but also tridiagonal; i.e., 



b ao 

a b\ (i\ 

ai b 2 



-J\f 



c 
-c ci 

-ci 



(0.25) 



In section 6 and 7, we show that the maps 0(_) and Or+) only involves three 
steps, in the following sense: 



—P2n(0,z) 
«2n 



I02n 
C2n 



( c 2 „ 2n \ 

-Can-lP2r»-l(0,z) H -Q^ & i)P2n(0, z) + C 2n p 2 n+1 (0, z) 

V fl 2» / 



{(3 = 1) (0.26) 



P2n(0, Z) 
P2n+l(0,z) 



-c 2 



— c 2r . 



n-lj ?2n-2(0) z ) + \/ a 2nC 2n ?2n^(0> Z ) 

V c 2n-2 



/ Q2n-2 ^(+) 
C2n-2 



2ri 



0-2n V ° 2n 

03 = 4) 



(0.27) 



The abstract map Or\ for t = appears already in the work of E. Brezin 
and H. Neuberger 0. This has been applied in jjL4j] to a problem in the 
theory of random matrices. 



1 Splitting theorems, as applied to the Toda 
and Pfaff Lattices 



In this section, we show how each of the equations 



drrioo { dmoo 
— — = A moo and — - — 
at; ati 



XI) 
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lead to commuting Hamiltonian vector fields related to a Lie algebra splitting. 
First recall the splitting theorem, due to Adler-Kostant-Symes [p]] and the 
R- version to Reiman and Semenov-Tian-Shansky |I5" |. The R- version allows 
for more general initial conditions. 

Proposition 1.1 Let g = k + n be a (vector space) direct sum of a Lie 
algebra g in terms of Lie subalgebras k and n, with g paired with itself via 
a non- degenerate a,d-invariant inner product (,); this in turn induces a 
decomposition g = k^ + n -1 and isomorphisms g ~ g* ; k -1 ~ n*, n -1 ~ k*. 
7Tk and 7r n are projective onto k and n respectively. Let Q, Qk and Q n be the 
groups associated with the Lie algebras g, k and n. Let 1(g) be the Ad* ~ 
Ad-invariant functions on g* ~ g. 

(i) Then, given an element 

eGg: [£,k] ck 1 and [e, n] C n , 

the functions 

v(t + Olic^ Wlth v e J (g) and i' e k x , (1.2) 

respectively Poisson commute for the respective Kostant-Kirillov symplectic 
structures of n* ~ k ± ; the associated Hamiltonian flows are expressed in 
terms of the Lax pairsf\ 

i = [-7r k V^(0, £] = fa^(0, £] for£ = s + e,Ce k 1 - (1.3) 

(ii) JTie splitting also leads to a second Lie algebra g R , derived from g, 
such that g* R ~ g R , namely: 

gR ■ [x,y]R = ^[Rx,y] + -z[x,Ry] = fax, 7r k y] - fax, -K^y] , (1.4) 
with R = 7r k — 7r n . The functions 

respectively Poisson commute for the respective Kostant-Kirillov symplectic 
structures of g R ~ wi/i i/ie same associated (Hamiltonian) Lax pairs 

e = [-vr k V^U] = KV^(£U] for £eg R . (1.5) 



6 (Ad s X ; y) = (X, Ad ff _iY}, .9 e G, and thus (\z,x],y) = (x,-[z,y\). 
7 Vip is defined as the element in g* such that d<p(£) = (Vip, d£), (gg. 
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Each of the equations (1.3) and (1.5) has the same solution expressible in 
two different ways^: 

£(t) = Ad K{t) £ = Ads-i( t) £o, (1.6) 

witlf\ 

K{t) = 7r gk e tv ^ o) , and S(t) = Ttg n e tv ^ \ 

dm 

Example 1: The standard Toda lattice and the equations — — = A l m 

.. dti 

for the Hankel matrix m^. Since, in particular, the matrix is sym- 
metric, the Borel decomposition into lower- times upper-triangular matrix 
must be done with the same lower-triangular matrix S: 

moo = S- 1 S T - 1 . (1.7) 

In turn, the matrix S defines a wave vector and operators^ L and M, the 
same as the ones defined in (0.22), 

1 oo 

y(t,z) '.= ett~ Uti Sx, L:=SAS-\ M := S(d+- j^A*" 1 )^ 1 , (1.8) 

2 i 

satisfying the following well-known equations^ 

d 

= = — with \L,M} = 1, 

oz 

dt n ~ 2 {L)b ° b dt n ~2 {L)s ^ 

K =1 2 [m *' L] wr\ [{Ln) ^ M] - (L9) 

The wave vector \t> can then be expressed in terms of a sequence of r-functions 
r n (t) = det m n (t), but also has the simple expression in terms of orthonormal 

8 naively written Ad K{t) £ a = K(t)£ K(t)-\ Ad s -i£o = S r_1 (i)^o-S , (t). 
9 with regard to the group factorization A = ng k A irg^A. 

10 In the formulas below x( z ) = ( z °i z ; z2 > ■ ■ ■) an d ^ ^ s the matrix such that j^x( z ) — 
dx{z)- 

n wliere the () s k and (){,<, refers to the skew-part and the lower-triangular (Borel) part 
respectively; i.e., projection onto k and n respectively. 
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polynomials, with respect to the moment matrix m,^: 

T n (t- [z- 1 ]) 



V(t,z) = e^E* 



z n 



r n (t)r n+1 (t) 

= e^^(p n (t,z)) n>0 . (1.10) 



The vector fields (1.9) on L are commuting Hamiltonian vector fields, in 
view of the Adler-Kostant-Symes splitting theorem (version (i)), 

8L trL i+1 

7 - = [-^VH l ,L] = [n n Vn i ,L}, H t = -—, VHi = L\ (1.11) 

Oli I + 1 

with 

L = A T a + b + aA, a and b diagonal matrices (1-12) 
for the splitting of the Lie algebra of semi-infinite matrices 

D — gloo — k + n := {skew-symmetric } + {lower-triangular} 

= k -1 + ir 1 := {symmetric } + {strictly upper-triangular}, 

(1.13) 

with the form (1.12) of L being preserved in time. Note that the solution 
(1.6) to (1.5) in the AKS theorem is nothing but the factorization of 
followed by the dressing up of A. 

Example 2: The Pfaff lattice and the equations — — = tCm + mA T . 

Oti 

Throughout this paper the Lie algebra T> = gl^ of semi-infinite matrices 

is viewed as composed of 2 x 2 blocks. It admits the natural decomposition 

into subalgebras: 

V = V-®V ®V+ = £>_ © £> ~ © £>+ © £>+ (1.14) 

where V has 2x2 blocks along the diagonal with zeroes everywhere else 
and where T> + (resp. T>_) is the subalgebra of upper-triangular (resp. lower- 
triangular) matrices with 2x2 zero matrices along V and zero below (resp. 
above). As pointed out in (1.14), V Q can further be decomposed into two Lie 
subalgebras: 

T>q = {all 2x2 blocks G T>q are proportional to Id} 

V+ = {all 2x2 blocks G V have trace }. (1.15) 
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Remember from (0.10) and (0.11) in the introduction, the matrix J and 
the associated Lie algebra order 2 involution J . The splitting into two Lie 
subalgebrasQ 

V = k 



n. 



1.16) 



with 
k = 



algebra of 



Q2n,2n 
Q2n,2n 



Q2n+2,2n+2 





Q 



2n+2,2n+2 



11 



V 

{a G V, such that J a = a} = {b + Jb, b G V} = sp(oo), 



• 7 

;i.i7) 



with corresponding Lie groups^ and Q n = Sp(oo), will play a crucial 
role here. Let and 7r n be the projections onto k and n. Notice that 
n =sp(oo) and Q n = Sp(oo) stand for the infinite rank afhne symplectic 
algebra and group; e.g. see |Tl[. Any element a G V decomposes uniquely 
into its projections onto k and n, as follows: 

a = iT^ci + 7r n a 

= | (a- - Ja+) + - (a - Ja )\ + |(a+ + Ja + ) + ^ (a + Ja )\ . 

(1.18) 

The following splitting, with 

k + = T> + + T?q~ and n + = n, 
will also be used in section 2; the projections take on the following form, 
a = vri c+ a + 7r n+ a 

= |(a+ - J a J) + ^ (a - Ja )^ + |(a_ + JaJ) + ^ (a + Ja )\ . 

(1.19) 



12 Note n is the fixed point set of J . 

13 (?k is the group of invertible elements in k, i.e., invertible lower-triangular matrices, 



with non-zero 2x2 blocks proportional to Id along the diagonal. 



13 



A v M : Pfaff /matrix integrals March 8, 1999 #1 §1, p. 14 



Note J intertwines and 7Tk + : 

J** = vr k+l 7. (1.20) 

For a skew-symmetric semi-infinite matrix m.oo, the skew-Borel decompo- 
sition 

moo := Q~VQ- 1T with Q E k , (1.21) 

is unique, as was shown in ||. Here we may assume moo to be bi- infinite, as 
long as the factorization (1.21) is unique, upon imposing a suitable normal- 
ization. Then we use Q to dress up A: 

L = QAQ- 1 . 

Then letting moo run according to the equations dm/dti = A l m + mA , 
we show in the next proposition and corollary that L evolves according to 
a system of commuting equations, which by virtue of the AKS theorem are 
Hamiltonian vector fields; for details, see 0. 

Proposition 1.2 For the matrices 

moo := Q 1 JQ~ 1T and L := QAQ^ 1 , with Q E Gk, 

the following three statements are equivalent 

(ii) V + ^-Q-' e n 

/•••i ■ . T i 

(in) — — = A m + mA . 

Whenever the vector fields onQ orm satisfy (i), (ii) or (Hi), then the matrix 
L = QAQ~ X is a solution of the AKS-Lax pair 



1 = 1-^^1 



7r n Z7, L 



Proof : Written out and using (1.18), proposition 1.2 amounts to showing 
the equivalence of the three formulas: 

(I) ^-Q- 1 + ((L% - J(L%) T J) + \ ((L% - J((V) ) T J) = 

14 
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din 

(III) A*m + mA T - — = 0. 

dti 

The point is to show that 



(D + = o. {i)- = (U)- = -J(n + )~ j. (Uu = |(n) , 



Q-^JQ- 1 = (III). (1.22) 
The reader will find the details of this proof in B . ■ 

2 Wave functions and their bilinear equations 
for the Pfaff Lattice 

Consider the commuting vector fields 

dm^/dti = A'moo + mooA 1 "* (2.0) 
on the skew-symmetric matrix m^t) and the skew-Borel decomposition 

rMt) = Q~\t)JQ T ~\t), Q(t) e £ k ; (2.1) 

remember from (1.17), Q(t) G means: Q(t) is lower-triangular, with along 
the "diagonal" 2x2 matrices c 2n 7, with c 2n 7^ 0. 

In this section, we give the properties of the wave vectors and their bilinear 
relations. Upon setting 

Qi = Q(t) and Q 2 = JQ T -\t), (2.2) 

the matrix Q(t) defines wave operators 

W 1 (t) = Q^eSr*^, W 2 (t) = Q 2 (t)e-^T^ Ti = JW{ 1T (t), (2.3) 

L-matrices 

L := L x := Q 1 AQ^ 1 , L 2 := -J(L t ) = Q 2 A T Q 2 ~ 1 , (2.4) 

15 



A v M : Pfaff /matrix integrals March 8, 1999 #1 §2, p. 16 



and wave and dual wave vectors 

z) = w 1 (t) x (z) *;(t, z) = wrWxiz- 1 ) = -.m 2 (t, z- 1 ) 
Mt,z) = w 2 (t) x (z) ** 2 (t,z) = w 2 -\ty x (z- 1 ) = JMt,^ 1 ). 

From the definition, it follows that the wave functions have the following 
asymptotics 

*i,2n(*, z) := e^ zk z 2n c 2n (t)^ 2n (t, z), ^ 1)2n = 1 + O^" 1 ) 

Vi^+ifaz) = e^ zk z 2n+1 c 2n (t)^ li2n+1 (t,z), Vl,2n+1 = 1 + 0(*" 2 ) 

*2,2n(*,*) = e^M-'^+l^l^^^ ^ = 1 + 

*2,2„ + l(t^)=e-5:^-^2n(-c^(t))^ 2)2n+1 (i,z), ^2,2n + l = l + 0(^ 2 ), 

(2.6) 

where the q are the elements of the diagonal part of Q. 
Theorem 2.1 The Qi, Li and satisfy the equations 

^ = -(7T k ^i)<5i ^ = -(J(^L[))Q 2 = (7r k+ 4)g 2 (2.7) 

^ = [-7r k Ll, Li] H = [7r k+ L 2 , L 2 ] (2.8) 
Li^i = L 2 # 2 = 2 _1 * 2 , (2.9) 

^ = ^ = -(L* - 7T k+ L 2 )* 2 = -(7T n+ L 2 )* 2 , (2.10) 

wift satisfying the following bilinear identity for all n,m E Z, 

c dz [' dz 

f * 1)n (i, z)#2, m (t', 2; -1 )- — — + f tt 2>n (f, z)V hm {t', z- 1 )—— = 0. (2.11) 
Joo Aniz Jo Zmz 

For later use, we shall also consider the "monic" wave functions, with the 
factors c 2n (t) removed, i.e., 

*i(t, z) := Q„ and * 2 (i, := Q ^ 2 (2.12) 

and the matrix L 1; normalized so as to have l's above the main diagonal, 
with Q := Q^Q, 
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U = Q L 2 Qo 1 = -QoJ(L 1 )Q 1 = -J(L 1 ) 



(2.13) 



Then, in terms of the elements of the matrix Q := Qq 1 Q, one easily 
computes by conjugation, that L\ has the following block structure: 



U = Qo'LiQo = (Qo 1 Q)HQo 1 Q)~ < 





Loi 








L\o 









* 


L>21 


L22 


L23 


* 


* 


L32 


L33 



with 



Q2i,2i-1 1 
Q2i+l,2i-l — <?2i+2,2i ~ <?2i+2,2i+l 





1 



Li 



* ~<?2i+2,2i+l _ <?2i+3,2i+l + <?2i+2,2i 



Theorem 2.2 Lj, Q, ^ 2 satisfy the following equations: 
^ = - (tf?)- - Qo 2 A(L-) + )\jQl) Q, 



and 
with 



_d_ 

dt n 



^Ht,z) = j(CLl) + + CL 1 l)o + Qo 2 Jm) + )Qt)^2it,z) 
= - ((L»)_ + (L5) + <&7((LJ)-)Qo 2 ) *2(«, *)■ 



(2.14) 



(2.15) 
(2.16) 
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The proof of Theorem 2.1 hinges on the following matrix version of the 
bilinear identities: 

Lemma 2.3 The matrices W\(t) and W 2 (t), defined in (2.3), satisfy 

W^W^t')- 1 = W 2 (t)W 2 (t')-\ (2.17) 
Proof : The solution to the equation (2.0) is given by 

m oo (t) = e^ t " Ak m oo (0)e^ t " ATk . 
Therefore skew-Borel decomposing m^it) and TOoo(O), we find 

q- 1 (o)jq t - 1 (o) = e -l>A i Q-i(t)jQT-i (t)6r i>A Tt (2.i8) 

and so, from the definition of W\ and W 2 , 

wr\o)w 2 (o) = g- 1 (o)jg T ^ 1 (o) 

= (Q(t)e^T^y 1 J (Q^eE^y- 1 using(2 .l 8 ) 

= w 1 (t)- 1 jw 1 (t) r - 1 

= W 1 (t)- 1 W 2 (t), (2.19) 

implying the independence in t of the right hand side of (2.19). Therefore, 
we have 

Wi(t) _1 W 2 (t) = W^t'^W^t'), for all t,t' G C°°, 

and so 

W^W^it') = W 2 {t)W^{t'). 

■ 

Proof of Theorem 2. 1 : The proof of equation (2.7) for Qi, namely 
follows at once from Proposition 1.2. 
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The proof of (2.7) for Q 2 = JQj 1 is based on the identity Jir^a 
n\L + Ja. Indeed, we compute 



—Q 2 = -JQ l —Q 1 Q 2 

= -JQj-'Qji^LlfQj-'Qjj 

= -J(ir k L\) T J 

= -J(ir k L{) 

= —Kk+JLl 

= ~^ k+ J(-JL 2 )\ using (2.4), 

= -n k+ j(-iy(jL 2 y 

Equations (2.8) and (2.10) for Li, L 2 and ^ 2 are then straightforward. 

Finally, the proof of the bilinear identity (2.11) proceeds as follows: By 
a well-known lemma (see [?]) , 



WiWimWitt')- 1 = I Vi(t,z) m(t',z) 

2 2 2 J oo 2 2 



dz 



2m z 

and so the statement of Lemma 2.3 yields 

£ * l( M) m\ z)± = j *,<«,*) 8 *5(f, 
whose (m, n)th component is 

Next we use the relations = -J^ 2 (t, z~ l ) and ^(t, = z ~ r ), 

to yield 

/" g?z f dz 

f *x(t, z) ® Jtf 2 (t', 2; -1 )- — — + i * 2 (i, z) <g> JV^t', z~ 1 )—— = 0, 
Joo Ztviz Jo Zmz 

which again componentwise leads to (2.11). ■ 
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Proof of Theorem 2.2 : To prove (2.15), remember from Theorem 2.1, 

^-Q- 1 = -vr k L" = -((L»)_ + J{LlYJ) - i((L») - J((^) ) T J); 

hence, taking the ( ) -part of this expression, yields 

dlogQo fdQ A 1 1 T 

Using the fact that Q , Qq 1 , Qo E Gk H X> commute among themselves and 
commute with J and the fact that T>qD± ) T>±Dq C X>±, we compute for 
Q = Qo X Q, U = Qo^iQo, (see (2.13)) 

^-Q- 1 = -Qo'QoQo'QQ^Qo + Q^QQ^Qo 



dt n 



— ~Qo l Qo + Qo l QQ 1 Qo 
= Qo'i-QoQo' + QQ'^Qo 

= Q^ 1 (-(L n 1 )_ + J(L n 1+ ) T j)Q 



= -{L n 1 )-+Qo 2 Jm) + yjQi 

Using this result and L 1 ^/ 1 (t, z) = z^>i(t, z), we find 

d$i(t,z) 
dt n 

= z n e ET^Q x ( z ) +e ^T^' (_(£«)_ + Q - 2 J((L n 1 ) + ) T JQl) Q x (z) 

= (l? - + g - 2 j((^) + ) T jg^i 1 (t, *) 

= ((L?) + + (L?)o + Qo 2 (^?) + )Oo)*i(*^)- (2-20) 

But, we also have that = Qq 1 ^^, z) and ^ 2 = Qo^2(t,z) satisfy, using 
W 2 = JW 1 ~ lT , 

= (Qo X Wr)- X {z) = {Qo 1 W 1 )\Q- Q 1 W 1 )- 1 i> 1 {t,z) (2.21) 
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9 ^ 2 d f ~ =(QoW 2 )- X (z) = (QoW 2 )iQ W 2 y\Q y 2 ) 

= (QoW 2 + Q W 2 )W 2 1 Q^ 1 (Q ^ 2 ) 

= (QoQo 1 + QoW 2 W^Q^)Q ^ 2 

= (QoQo 1 + QoJ(W 1 W 1 - 1 )Q 1 )Q ^ 2 

= (QoQ 1 + Q J(W 1 Wr 1 ) T JQo 1 )Qo^2 

= (-JQoQo 1 + J(Qo 1 WiW^ 1 Q ) T )JQ ^ 2 

= J(-Qo 1 Qo + Qo 1 WiW 1 - 1 Qo) T JQo^2 

= J \{Q^W l )\Q^W l )- 1 ){Q^ 2 ). (2.22) 

Comparing (2.21), (2.22) and (2.23), and using 

-3{LX) = LI 

and so, in particular, 

-jm)±) = (A% and - J((L?) ) = (L5) , 

9 -^§^ = "(A)- + (^)o + Ql( J{L n 2 )_)Q- Q 2 )%(t, z), 
which establishes theorem 2.2. ■ 



3 Existence of the PfafF r-function 

The point of this section is to show that the solution of the Pfaff Lattice can 
be expressed in terms of a sequence of functions r, which are not r-functions 
in the usual sense, but enjoys a different set of bilinear identities and partial 
differential equations. 

Proposition 3.1 There exists functions T 2n (t) such that 

V>i,2n{t,z) = 7- and ip 2 , 2n [t,z) = — . (3.1) 

T 2n {t) r 2n+2 (t) 

The proof of proposition 3.1 will be postponed until later. For future use, 
we define the diagonal matrix 

h = dicLg(...,h- 2 ,h- 2 ,h ,h ,h 2 ,h 2 ,...)eV , with h 2n = ^±i. (3.2) 

r 2n 



21 



A v M : Pfaff /matrix integrals March 8, 1999 #1 §3, p. 22 



Theorem 3.2 

*l,2n(M) 
*l,2n+l(*»2) 

*2,2n+l(*»*) 



3 ^i, 2 ' 2n/2n(< ~ *]) 



= e 



= e 



y / T2n(t)r 2 „+2(t) 
^ z 2n (z + d/dtjT^it-jz- 1 ]) 
\J T2n(t)T2n+2(f) 
- Uz-* z 2n+l T 2n+2(t + \z\) 
\J T2n(t)T2n+2(t) 

• 5>*-' 2n+i(£ -1 - <9/9ti)r 2n+2 (t + [*]) 



\l T2n(t)T 2 n+2(t) 



(3.3) 



mt/i £/ie r 2n (t) satisfying the following bilinear identity 

I T 2n (t - [z- l ])T 2m+2 (t> + [,-l])eE(*.-^'^-^ 

+ / r 2n+2 (t + [*])r 2m (t' - [^]) e E(^-* i )-- i ^n-2 ? n|i = (3 4) 
T/ien L /ias #ie following representation in terms of the Pfaffian r-functions: 



with (• = Wi ) 

( -(logr 2n )- 



n.n • 



Loo 










L±o 


L ^ 


L\ 2 





* 


L 2 i 


L 22 


-^23 


* 


* 


L32 


-^33 



T2n T2n + 2 



(logr 2n+2 )7 





1 



-<n+l,n 



* (logr 2n+2 )- 



(3.5) 
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The following bilinear relations are due to ||: 

Corollary 3.3 The functions r 2n (t) satisfy the following "differential Fay 
identity'^ 

{T2n(t - [u]),T 2n (t - [v])} 

+ ~ Ofa^ - [u])T 2n {t - H) - T 2n {t)r 2n {t - [u] - H)) 

= uv(u - v)r 2n _ 2 (t - [u] - [v])r 2n+2 (t), (3.6) 
and Hirota type bilinear equations, always involving nearest neighbours: 
( ~ 1 d 2 \ 

p k +i{d) - -wtk: T 2n o r 2n = p k (d) r 2n+2 o r 2n _ 2 (3.7) 



2 dt x dt k 



+3, 



k,n = 0,1,2, 



Lemma 3.4 Consider an arbitrary function ip(t, z) depending on t e C°°, 
z G C, having the asymptotics <f(t,z) = 1 + O(-) for z/oo and satisfying 
the functional relation 

V(t-[Z 2 \ Z1 ) = <p(t- \Zi\z 2 ) teC oo 
(f(t,Zi) V(t,Z2) 

Then there exists a function r(t) such that 

Proof : See appendix (Section 10) 

Lemma 3.5 The following holds for the Pfaffian wave function and ^ 2 , 
as in (2.6), 

V>l,2n(*- [Z2 l lzi) = ^l,2n(* ~ fa 1 ], Z 2 ) 
l/>l,2n(t, Zi) 1pl,2n{t,Z 2 ) 

and 

^2,2n- 2 (t - [z-\ Z- X )iJ h2n (t, Z) = 1. (3.9) 



l {f,9} = f'9-fg', where ' = d/dt x . 
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Proof : Setting (2.6) in the bilinear equation (2.11), with n \— > 2n, m 
2n — 2, yields 

C2n-2\t) Joo Zm 

Setting 

t-f = [^r 1 ] + [z, 1 ] 

in the above and using e^i — 1/(1 — x) yields 

C2n f ^l,2n(t,z)iJ 2 ,2n-2(t' ' , Z' 1 ) dz 



c 2n -2 Joe (i - Jl) (l - a) 2ni 



c 2n Jz=o V Z\) V z 2 J ' ' 2m 

Since the integrand on the right hand side is holomorphic, it suffices to eval- 
uate the integral on the left hand side, which can be viewed as an integral 
along a contour encompassing oo and the points z\ and z 2 , thus leading to 

^l,2n(t, ^l)^2,2n-2 (*~ [z^ 1 ] ~ [z^ 1 ] , Zf 1 ) = ^l,2n (*, Z 2 ) ^2,2n-2 (*~ [^f 1 ] ~ [zj 1 ] , 

(3.10) 

with 

^l,2n(*, *) = 1 + O (z- 1 ) , ^ 2 ,2n- 2 (* " [z?\ ~ [z 2 \ Z' 1 ) = 1 + 0(z' 2 ). 

Therefore, letting z 2 / oo, one finds 

1pl,2n(t, Zi)^ 2 ,2n-2{t ~ [z±\ Z^ 1 ) = 1, (3.11) 

yielding (3.9), and so, upon shifting 1 1— > t — [z 2 l ], 

1p2,2n-2(t ~ [Z^ 1 ] - ta" 1 ],^" 1 ) = ( . \ _u V 5 

similarly 

^2,2n- 2 (t - ft -1 ] - teW) = ,, \ -x, r - (3.12) 

Setting the two expressions (3.12) in (3.10) yields 

^l,2n(* ~ fe" 1 ], *l) _ ^l,2n(t - [z^], Z 2 ) 



^l,2n(Ml) ^l,2n{t,Z 2 ) 
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Proof of Proposition 3.1 : From Lemmas 3.4 and 3.5, there exists, for 
each 2n, a function r 2n such that the first relation of (3.1) is satisfied, i.e., 

T2n(t-[Z- 1 }) 



and so from (3.9) 

^2,2n-2(t - \z~\z~ 1 ) 

thus leading to 



T2n{t) 



1 



T2n(t) 



l/>l,2n(t,z) T 2n {t -[Z 
T2n(t+[z\) 



^2,2n-2(*, Z) = , 
r2n{t) 

which is the second relation of (3.1). ■ 

Proof of Theorem 3.2 : At first, remembering that Q = Q^Q, observe 
that 

e^ UZ \(Q)x(z))2n = (Qo ll MM))2n 

= e£^V> 1)2n (M) 

= c Y.Uz* z 2n T 2n(t - [Z- 1 ]) 
T2n{t) 

V ^! r 2n(t) J 

showing that a few subdiagonals of the matrix Q are given by 
/ • 



Q 



1 
1 




<?2n,2n-2 <?2n,2n-l 
92n+l,2n-2 92n+l,2ra-l 


1 
1 



with 



9 A S 2 (-d)r 2n 

Q2n,2n-1 ~ ~ T7~ >-OgT 2n , q 2n ,2n-2 ~ • 

oh T 2n 



(3.13) 
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Remembering that 



If>l,2 n = l + 0(Z- 1 ) 



*i, 2n+ i(t, z) = eI>*V n+ Vi,2n + i(t, z), Vi,2n+i = 1 + 0(z- 2 ) 

*2,2n(*,*) = e-5>*~V" +1 V> 2>2n (M), ^2,2n = 1 + O(z) 

*2,2n+l(*^) = -e-^ z ' k Z 2n ^ 2>2n+1 (t,z), ^ 2 , 2n+1 = 1 +0(Z 2 ), 



we now compute, using theorem 2.2, 



(3.14) 



(3.15) 



= (((L 1 ) + + (L 1 ) + Q 2 J(L 1+ ) T JQ 2 )y 1 (t,z)) 2n (3.16) 



In 



and 



z 2n+l ^ 2 , 2n {^z) 



d 



= {{j{{L 1 ) + + {L l ) + Qfj{L l+ yjQl))i> 2 {t,z)) 2n . (3.17) 
In this expression, the matrix equals, according to (2.13), 

(L 1 ) + + (L 1 )o + Q 2 J(L 1+ ) t JQ2 = 



30,-1 


1 














3i,-i - 320 


-321 


1 

















321 


1 








C o/ C 2 





331 - 342 


-343 


1 

















343 


1 








C l/ C 4 





353 — 364 


-365 ••• 



V 



and, acting with J on this matrix, 
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J ((L0+ + (LOo + Q 2 J(L 1+ ) T JQl) = 

( ■ \ 



-521 


1 














5i,-i - 520 


5o,-i 


C o/ C 2 

















-543 


1 








1 





531 — 542 


521 


C2/C4 

















-565 


1 








1 





553 — 564 


543 ••• 



V 



using the fact that 



J 

















0) 



Therefore the 2nth rows of both matrices respectively have the form 

(0,...,0,52n,2n-l(*), 1,0,0,...) 

T 

2n 

(0, . . . , 0, -52n+2,2n+l(t), 1, 0, 0, . . .) 

T 

2n 

and thus from (3.16) and (3.17), and the expansions (3.14) and (3.15), we 
have 

{jjty + *) ^' 2n ^' ^ = 22n,2n-l Wl,2n + #l,2n+l 
{jUj- ~ 4>2,2n(t, Z) = q 2 n+2,2n+l{t)lp2,2n + Z~ l ^ 2 ,2n+1 (3.18) 

and so, using the expression (3.13) for 52n,2n-i an d the first expression (3.1), 

Z 2n+1 ^,2n + l(t,z) 

= [z + ^ Z 2n ip li2n (t, z) - 52n,2n-l(t)^ n ^l,2n(t, z) 

= (z + z 2n ip 1M (t,z) + ( Ai og7 - 2n ( t )j z 2n if; 1M (t,z) 
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/ , S \ „ 2 „ T 2 „(t -[;-■]) id \ „ 2 „ T 2 „(t-[z-']) 

= [ Z+ Sh) Z r 2n (t) + [W^j Z r U t) 
and similarly, using the relation (3.17), 

2 „ 2w+1 (-^ 1 + lr) 7 "^( t + M) r , 9n x 

2 1p2,2n+l{t,Z) = Z — . (d.iUJ 

Therefore, we also have 

1 / <9 2 \ 



thus 



n 1 fez- a\ ° 2 \ -S2(d)r 2n 

?2n+l,2n — U, ?2n+l,2n-l — «J2i, — °V ~~ "^TJ T 2n — • ^O.^i J 

T2n \ Oh ) T 2n 

Setting n h-> 2n and m i— > 2n in the bilinear relation (2.11) and substitut- 
ing, using (2.6) and the expressions for ipi ;2 n(t, z ) an d i>2,2n(t, z) m the proof 
of proposition 3.1, 

*i,2n(c, = 2; c 2n (t) r— 

T 2 n{t) 

and 
into 

^ i / iTr / a „\,t, (4.1 -i\ dz 



f dz f dz 



yields 



C 2n(t) /■ c y-fa-t0^ r 2n(t- [g *] )r 2w+2 (g + [g X ]) 



C2n(i') • '°o T 2n (t)T 2n+2 (t>) 2mz 2 

28 



A v M : Pfaff /matrix integrals March 8, 1999 #1 §3, p. 29 



| C2 n (t) f c T.(^-t k )z- k T 2n+2(t+ [z])T 2n (f - [z]) dz 
C2n(t) JO T 2n+2 (t)T 2n (t') 2m' 

Setting t' — t + [a] amounts to replacing the exponential: 

e £<fc-f fc )** = l-az, eEtt-fc)*-* = — L_ 

1 — a/z 

so that the first integral has a simple pole at z = oo and the second integral 
one at z = a. Evaluating the integrals yield 

—a -7T7T + at = 

r 2 (+/\ T2n+2(t') ' 

c 2nV ) r2n(t /) 



i.e., 



T2n(0 

yielding the following relation, which involves a constant c n , independent of 
time, 

cl(t) = c n ^- = c n .h-^t). (3.22) 

Rescaling r 2n h- > r 2ri /(cic 2 • • • c n _i) yields (3.3). Using the expressions for 
ipi,2n(t,z) and ^2,2n(t,z) (see the proof of proposition 3.1), (3.19), (3.20), 
(3.21), (2.6) and substituting (3.3) into (2.11) yields (3.4). 

Finally to derive the form (3.5) of the matrix L, set (3.13) and (3.20) in the 
elements just below the main diagonal of matrix (2.14), to yield (• = d/dti) 

— ?2n,2n-l ~~ 92n+l,2n-l + <?2n,2n-2 

T2n\ 2 _ (^(-9) - + S 2 {-~d)r 2n 



- T 2n ' T 2n T 2r, 

T2n [T 2 r. 



and 



T~2n \T2r, 

= (logr 2n )- 



92n+l,2n-l — Q2n+2,2n — 



(S 2 (-d) - §z)T 2n S 2 (-d)r 2n+2 

T 2n T 2n+2 

_S 2 {dyr 2 n _ S 2 {-d)T 2n+2 

j 

T~2n T 2n+2 
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concluding the proof of theorem 3.2, upon substituting these relations into 
(2.14). ■ 



4 Semi-infinite matrices moo, (skew-) orthogonal 
polynomials and matrix integrals 

4.1 dm/dtk = A fc m, orthogonal polynomials and Her- 
mit ean matrix integrals. 

For the sake of completeness and analogy, we add this subsection, which 
summarizes some ofH. Consider a t-dependent weig ht Pt (dz) := e~ Vt ^dz : = 
e - v ( z )+J2 t i zl dz = e£^* z *p(dz) on R, as in (0.0) and the induced t-dependent 
measure 

e Tr(-V(X)+5>*) dX> 

on the ensemble 7i n of Hermitean matrices, with Haar measure dX; the latter 
can be decomposed into a spectral part (radial part) and an angular part: 

n 

dX := J] dXu ]J (dftXij dSsXij) = A 2 (z)d Zl ■■■dz n dU, (4.2) 

1 ^<i<j<n 

where A(z) = (zi — Zj) is the Vandermonde determinant. Here we 

X<i<j<n 

form the following matrix integral 

f ( ?H-viX) + j:u*) dx = Cn [ A \z)f[ Pt (dz k ). (4.3) 

The weight Pt(dz) defines a (symmetric) inner product 

(f,9) sy = f f(z)g(z)p t (dz) 



and so, the moments 

fJUfi) := (z\z j ) s y= I z k+t e^ zi p(dz)=p i+£d (t) 

satisfy 

_ f z iW e E^ k p (dz). 

JR. 



dt. 
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Therefore the semi-infinite moment matrix moo = {Hij)i,j>o satisfies 

^ = AV^ = rrwV^ (4.4) 
dti 

The point now is that the following integral can be expressed as a determinant 
of moments, namely 



(*)=/ e-WWdX = I A\z)Hp t (dz k ) 

,. n 

= / R .i:<iet(^ )Z -i) i<u<i> np«(^) 

aGSn k = l 

= J s det (4w~ 2 )i <w< „ n pt(dz„{t)) 

<j£S n k=l 



aes n 

n\det([ z e+k ' 2 p t (dz)) 

V-/R / l<£,k<n 

n\ det(// i j)o< i j< ri _i 



l<e,k<n 



is a r-function for the KP-equation; also in view of (4.4) and the upper-lower 
Borel decomposition (0.3) of moo, the integrals form a vector of r-functions 
for the Toda lattice. 

4.2 dm/dtk = A fc m+mA Tfc , skew-orthogonal polynomials 
and symmetric and symplectic matrix integrals. 

Consider a skew-symmetric semi-infinite matrix 

ra<x>(0 = (A*«(*))tj>o, with m n (t) = (£))o<ij<n-i, 

satisfying 

dm^/dtk = A k moo + mooA Tn . (4.5) 

Then we have shown in previous sections that, upon skew-Borel decompos- 
ing moo, these equations ultimately imply the existence of functions r(t) 
satisfying the bilinear equations (3.4). Remember also 

h{t) = diag(. .. ,h- 2 ,h-2,h , h ,h 2 ,h 2 , . . .) G V , with h 2n (t) = T2n+2 ^ , 

7"2n(lJ 
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Here, we need the Pfaffian pf(A) of a skew-symmetric matrix A = (aij)o<i,j<n-i 
for0 even n: 

pf(A)dx A ■ • • A dx n _i 
1 



n! 



2 n/2 



ctijdxi A dxj 

0<i<j<n—l j 
' n /2)\ £ ( a ) a io,ii a i2 



i 3 ■ • • o>i n _ 2t i n _ 1 J dxo A ■ ■ • A dx n -i, 

(4.6) 



so that pf(A) 2 = det A. We now state the following theorem due to Adler- 
Horozov-van Moerbeke^. 

Theorem 4.1 Consider a semi-infinite skew- symmetric matrix moo, evolv- 
ing according to (4-5); setting 

T2n(t) =pf(m 2n (t)), and h 2n = P ^™ 2n+2 ®\ (4.7) 

pf{m 2 n{t)) 

then the wave vector defined by (3.3) is a sequence of polynomials, except 
for the exponential, 

*i, k (t,z) =e^ Uzl q k (t,z), (4.8) 
where the q k 's are skew- orthonormal polynomials of the form (0.17) satisfying 

((Qi, Qj) sk )o<i,j<oo = J, with (y\ z j ) sk := fuj. (4.9) 

The matrix Q defined by q(z) = Qx( z ) i> s the unique solution (modulo signs) 
to the skew-Borel decomposition of m^: 

m OB {t) = Q- 1 JQ T - 1 , withQek. (4.10) 

The matrix L = QAQ^ 1 , also defined by 

zq(t,z) = Lq(t,z), 

and the diagonal matrix h satisfy the equations 

— = -vr k L\ L . and h' 1 — = 2n k (U) . (4.11) 
oh 1 1 Otj 



15 In the formula below (zq, z'i, . . . , i n -2> hi-i) = c(0, 1, . . . , n — 1), where a is a permu- 
tation and e(cr) its parity. 
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Sketch of proof, at first note that looking for skew-orthogonal polyno- 
mials is tantamount to the skew-Borel decomposition of m^, so that (4.9) 
and (4.10) are equivalent. The skew-orthogonality of the polynomials (0.17) 
follows from expanding the determinants explicitly in terms of ^-columns, 
upon using the expression for the pfaffian in terms of a column 

£ (-l) k a kiP f(0, k, ...,£- 1) = p/(0, 1, i). 

o<k<e-i 

For details, see ||. On the other hand, Theorem 3.2 gives ty(t, z) and hence 
Q in terms of r n (t) of (4.7). By the uniqueness of the decomposition (4.10), 
the two ways of arriving at Q, (0.16) and (3,3) must coincide. ■ 



Important remark: The polynomials (0.16) provide an explicit algorithm to 
perform the skew-Borel decomposition of the skew-symmetric matrix moo. 
Namely, the coefficients of the polynomials qt provide the entries of the matrix 
Q. This fact will be used later in the examples. 



Symmetric matrix integrals Here we shall focus on integrals of the type 

/ e^C-^W+Er^dX, (4.12) 
where dX denotes Haar measure 

dX := J] d ^ X H = 1^)1^1 ■■■ dz n dU, (4.13) 

l<i<j'<n 

over the space S2 n of symmetric matrices. As will appear below, the integral 
(4.12) leads to a skew-inner-product with weig ht p t {z)dz := e- v *&dz : = 
e - v ( z )+Y, tiZl dz = e^'* 2 ' p{z)dz on an interval C R, as in (0.1), 

(f(x),g(y)) := J f(x)g(y)e(x - y)p t (dx)p t (dy) (4.14) 

and therefore skew-symmetric moments 16 



Pij{t) = 




x % y 3 e(x - y)pt(x)p t (y)dxdy 



e{x) = 1, for x > and = — 1, for x < 0. 
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(x l y j - x 3 y l )p t {x)p t {y)dxdy 

x>y 

(F j (x)G i (x) - Fi{x)GAx)) dx. (4.15) 

R 

where 

/X 
y l p t {y)dy and G^x) := F-(x) = x l p t (x). 
-oo 

By simple inspection, the moments pkiif) satisfy 
dp k e 



0f j /r2 (x k+ Y + x k y l+l )e{x - y)p t (x)p t (y)dxdy 

— Hk+i,e + Pk,e+i, 

and so satisfies (4.5). 

According to Mehta |TB[, the symmetric matrix integral can now be ex- 
pressed in terms of the pfaffian, as follows, taking into account a constant 
c 2n , coming from integrating the orthogonal group: 



1 / e Tr (-v(x) + E^) dX 
2n)\ Jso n (E) 



(2n)\ Js 2n (E) 

|A 2n (^)| Y[pt(z,i)dz, 



i=l 

2n 



(2n)\ Jn 2n 

/ det \z\ , -iPtizi+i) ) W dzi 

n 

II Pt{z 2k )dz 2k 



CO<Z2<Z4<---<Z2n<00 



k=l 



Z'Z . . fZ2n 



det / z\p t {z^dz x 

j Z 2 ) • • • ; / Z 2n— \Pt ( z 2n— l)dz 2n _ 1 > ^2n 



oc 

ft 



0<i<2n-l 



II Pt(z2k)dz 2 k 

CO<Z2<Z4<---<Z2n<00 



det (Fj^) , 4 , -^(^4) - Fi(z 2 ) , z\ , . . . , Fj(02n) - ^(^-2) , 4 

n 

IJcfe; det^F;^) , Gi(z 2 ) , i^n) , G,- t {z 2 



0<i<2n-l 



-00<Z2<24<---<^2n<00 



'0<i<2n-l' 
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1 f n 

= ^|7 Rn n^ de*>( F i(yi) , G i(yi) , F i(yn) , ^ ''<(//). 

VR / 0<i,j<2n-l 

using de Bruijn's Lemma ||13||,p.446, 

= Pf (/ i ij)o<i,j<2n-l 

= r 2n (t), 
which is a Pfaffian r-function. 

Symplectic matrix integrals: Here we shall concentrate on integrals of 
the type 

l e 2Tr (-V(X) + £^) rfX; ( 4.1 6 ) 

where dX denotes Haar measure^] 

N - 

dX = Y[ dX k J[ dX M dX kl dX kl dX kl , 

1 k<l 

on the space T 2 n of self-dual N x N Hermitean matrices, with quaternionic 
entries; the latter can be realized as the space of 2N x 2N matrices with 
entries xfj £ C 



'(0) v (i) 



T 2 n — I X — (Xke)i<k,e<N, Xke — ( My -\o) ) with — Xj e , 

A more exotic skew-symmetric matrix satisfying (4.5) is given by the 
moments, with V(y,t) = e 2 ^ v ^ + ^ taya \ 

IMS® = I {y\y^- y ^ + ^y a )l E (y)dy 

= I {y>e v ^\yie v ^}I E (y)dy 
Jr 

= f (Gi(y)Fj(y) — F i {y)Gj{y))dy 1 (4.17) 

JR. 

upon setting 

Fj{x) = x ] e v{x ' t) and Gj{x) := Fj(x) = (a? e v . 

17 x means the usual complex conjugate. The condition on the 2x2 matrices Xke 
implies that Xkk = Xkl, with X^ € R. and / the identity. 
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That moo satisfies (4.5) follows at once from the first expression (4.17) above. 



m(t) = J {y k ,y e }p t (y) 2 dy 

= J{k-i)y k+l - l p t {yfdy 
^ = 2j{y\yt}y^y) + E^) dy 

= J((k + i- £)y k+i +^ + (k — £ — Oy*^" 1 ) Pt(y) 2 dy 

= H>k+i,i + pk,l+h 

thus leading to (4.5). Using the relation 

n - s,-) 4 = det (x\y 4 (4)'... 4 Kr) o<t<2ni , 



l<i,j<n 



one computes, using again de Bruijn's Lemma, 

J_ f e 2Tr(-V(X)+^U^) dX 

{n)\ Jr 2 „ 



1 r n 

= n (^^rn^"'""* 

U. JK i<jj< n i=1 
1 r n 

det (4 (*•)' 4 (4)' • • • 4 (4)') 

i /■ n 

= -. Ud yi dettaj/i) G t ( yi ) ... F t (y n ) G^y 
= det 1 / 2 (JjGiMFjiy) - F i (y)G j (y))dy 



0<i<2n-l 



0<i<2n-l' 



0<i,j<2n-l 



j<2n-l 

= T 2n ( t ), 

which is a Pfaffian r-function as well. 
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5 A map from the Toda to the Pfaff lattice 

Remember from (0.1), the notations pt(z) = p(z)e^ tkzk , and p'/p = —g/f. 
Assuming, in addition, f(z)p(z) vanishes at the endpoints of the interval un- 
der consideration (which could be finite, infinite or semi-infinite), one checks 
the the t-dependent operator in z, 

17 d r— 
n t ■ = \ — rvfPt 



Pt dz 



= ^/W-^-pW. ™thg t (z)=g(z)-f(z)Y:kt k z k -\ 

(5.1) 

maintains H + = {1, z, z 2 , ...} and is skew-symmetric with respect to the t- 
dependent inner-product ( , ) s t y , defined by the weight p t (z)dz, 

(n t (p,?p) s t y = / (n t (p)(z)ip(z)p t (z)dz = - (p(n t ip)p t dz = -{(p,n t ip) s t y . 
Je Je 

The orthonormality of the t-dependent polynomials p n (t, z) in z imply 

(Pn(t,z),p m (t,z)) S t y = 5 mn . 

The matrices L and M are defined by 

zp = Lp and e ~^ tkzk -^-e^ tkzk p = Mp. 

dz 

The skewness of n t implies the skew-symmetry of the matrix 

j\f(t) = f(L)M - t^-?-(L), such that n t p(t, z) = J\fp(t, z); (5.2) 

so j\f(t) can be viewed as the operator n t , expressed in the polynomial basis 
(Po(t,z),p!(t,z),...). 

In the next theorem, we shall consider functions F of two (non-commutative) 
variables z and n t so that the (pseudo)-differential operator in z and the ma- 
trix 

u t :=F(z,n t ) and U := F(L,M), 
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related byQ 

F(z,n t )p(t,z) = F(L,X)p(t,z), 
are skew-symmetric as well. Examples of F's are0 

F(z,n t ) :=n t , n^ 1 or {/,nf +1 }t, 

corresponding to 

F{L,N) = M, AT" 1 or {Af 2k+ \L e y. 

Theorem 5.1 Any Hdnkel matrix evolving according to the vector fields 

drriooit) k 

— — = A moo 

ot k 

leads to matrices L and M, evolving according to the Toda lattice (1.9). 
Consider a function F of two variables, such that the operator u t := F(z, n t ) 
is skew- symmetric with respect to ( , ) s t v and so the matrix 

U{t) = F(L(t),Af(t)), defined by u t p(t, z) =Up(t,z) 

is skew- symmetric. This induces a natural lower-triangular matrix 0(t), 
mapping the Toda lattice into the Pfaff lattice: 



p n (t,z) = {S(t)x(z)) n orthonormal with respect to 



Toda lattice < 



0<i,j<oo 



S^S 



1 cT-1 



L(t) = SAS 1 satisfies 



dL 



L j ,L 



J = 1,2, 



18 with the understanding that F(L,Af) reverses the order of z,u in F(z,u) 
19 {A, B}^ = AB + BA. 
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map 0{2t) such that < 



-U(2t) = 0^ 1 (2t)JO T ' l {2t) 
0{2i) is lower-triangular 

[ 0(2t)S(2t) e g k 



q n (t, z) = (0(2t)p(2t, z)) n , skew-orthonormal with regard to 
m^t) := -S- 1 (2t)U{2t)S T - 1 {2t) = Q' 1 (t)JQ T ' 1 (t) 



Pfaff lattice < 



0<i,j<oo 



dL 

L(t) := 0{2t)L(2t)0(2t)- 1 satisfies — = [-ir k L j ,L] , j = 1, .. 

3 

Proof : Since U(t) is skew-symmetric, it admits a skew-Borel decomposi- 
tion 

-U{t) = 0~\t)JO r -\t), with lower-triangular 0(t). (5.3) 
But the new matrix, defined by 

m^it) := -S- 1 {2t)U{2t)S T - 1 (2t), (5.4) 

is skew-symmetric and thus admits a unique skew-Borel decomposition 

rhoo(t) = Q- l (t)JQ(ty- 1 with Q(t) e (5.5) 

Comparing (5.3), (5.4) and (5.5) leads to a unique choice of matrix 0(t), 
skew-Borel decomposing —U(2t), as in (5.3), such that 



0(2t)S(2t) = Q(t) e g k . 



(5.6) 
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Using 



and 



we compute 

drhoo 



Otk 

^(2t) = -(n bo L k (2t))S(2t) : 



dt 



k 



it) 



S- 1 ^-{2t)S- 1 U{2t)S T -\2t) - S-\2t) L^-U{2t)\ S T ~ 1 (2t) 



i-^dS , ,n. c ,t~i 



+S- 1 (2t)U(2t)S T - 1 ^— (2t)S 
ot k 

= -S- 1 (ir bo L k (2t))US T - 1 - S- 1 [ir sy L k ,U]S T - 1 - S- l U(ii bo L k ) T S T ' 1 
= -S-\ix bo L k + n sy L k )US T - 1 - S- l U{{u bo L k y - HsyL^S^- 1 
= -S- 1 L k US T - 1 - S- 1 UL Tk S T -\ using (5.6) below 
= -A h S- 1 US T - 1 - S- 1 US T - 1 A Tk S T S T - 1 , using L k = SA k S~\ 
= A k rh 00 (t)+m 00 (t)A Tk . 

For an arbitrary matrix A, we have 

A = A T A = (A bo ) T - A sy . (5.7) 

Indeed, remembering thatQA^ = 2A_ + Aq and A sy = A + — A_, one checks 

(A bo ) T -A sy -A = 2(A_) T +A -(A + -A_)-A_-A + -A = -2(A + -(A_) T ). 

so that the left hand side vanishes, if the right hand side does; the latter 
means A is symmetric. 

We now define L{t) by conjugation of L{2t) by 0{2t): 

L{t) := 0(2t)L(2t)0(2t)- 1 = 0{2t)S{2t)AS-\2t)0{2t)- 1 = Q(t)kQ-\t). 

Therefore the sequence of polynomials 

q(t,z) := Q(2t)p(2t,z) = 0(2t)S(2t) X (z) = Q(t) X (z) 

20 A± means the usual strictly uppcr(lower)-triangular part and Ao the diagonal part in 
the common sense. 
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is skew-orthonormal 



(qi(t,z), qj (t,z)) sk = J i:i 



with regard to the skew inner-product specified by the matrix m^: 

(z\zi)f = /2y(f). 
In the last step, we show that ((p,?p) sk = (</?, \wp) sy . Since 

U(2t) = -0' 1 (2t)JO T - 1 (2t), (5.8) 

we compute 

( qi (t,z),(u 2tq ) 3 (t,z)) s £ = ((Op) i (2i),(uOp) J (20)^ 

= ((Pp) i (2t),(Oup) i (20>^ 
= <(Pp) i (2t),(OWp) J -(20>^ 
= (0(2t) (^(2t),p £ (2t))^>o(OW) T (2t)) 4J 

= (0(2t)/(OW) T (2t)) y . 

= (0(2t)U T (2t)0 T (2t)) i:j 

= -(0(2t)W(2t)0 T (2t)) 4J 

= using (5.8). (5.9) 

Therefore, defining a new skew inner-product ( , ) sk ' 
we have shown 

and so by completeness of the basis qt, we have 

( , )f = ( , )f , 

thus ending the proof of Theorem 4.1. ■ 
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6 Example 1: From Hermitean to symmetric 
matrix integrals 

Striking examples are given by using the map 0(t) obtained from skew-borel 
decomposing A/" _1 (t) and N{t); see (5.2). This section deals with M~ x {t), 
whereas the next will deal with M{t). 

Proposition 6.1 The special transformation 

U{t) = M-\t) = (f(L)M - (2f) 

maps the Toda lattice t -functions with initial weight p = e~ v , V = —g/f 
(Hermitean matrix integral) to the Pfaff lattice r-functions (symmetric ma- 
trix integral), with initial weight 



Pt\z) :-- 



To be precise: 



-±(V(z)+\o g f(z)~2j2?Uz') = . e -V(z)+J2?Uz 



Toda lattice 



p n (t, z) orthonormal polynomials in z for the inner-product 

(<P,ip}t V == J V{z)^{z)pt{z)dz, 
fiij(t) = (z\z 3 ) s t y and m n = (fHj)o<i,j<n-i, 
r„(t) = det m n = 1 / e Tr ^ x ^T **)dX 



-M-\2t) = 0-\2t)JO T -\2t) 
map 0(2t) such that ^ 0{2t) is lower-triangular 

0(2t)S(2t) G £ k 



21 Remember p'/p = —V = —g/f. 
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q n (t,z) — 0(2t)p n (2t, z) skew- orthonormal polynomials 

in z for the skew-inner-product (weight p), 



Pfaff lattice < 



2t 



1 



<p(x)ij)(y)e(x - y)pt{x)p t (y)dx dy 



2 J Jr2 

p,ij(t) = (x\y J )f andm n = (/%)o<ij<n-i 



f 2n (t) = pf{m 2n ) = J- / e^W+Er^dX. 

In the first integral defining T n (t), dX denotes Haar measure on Hermitean 
matrices (see section 4.1), whereas the second integral f 2n (t) involves Haar 
measure on symmetric matrices (see section 4.2) 

Proof : At first, check that 



Indeed, 



did 
dx 1 dx 



/d 
S(x - y)(p{y)dy using — e(x) = 25 (x) 



ip(x). 



Consider now the operator 



— TVfPt) ' so that ™ t p = n^p = Af^p, 



u t = n t = 

according to (5.2). Let it act on a function <p(x): 



n t <p{x) = 



1 



d 



-i 



y/f(x)p t (x) \dxj \ f(x 

i e(x - y) Pt{y) 

* sjf{x)p t {x) 2 \ f(y) 
43 
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One computes 



(<p>if>) 



(ip, n 2t ip) s 2 y t 



{<P, n 2tV>2i 



2 J JR.* 
1 



\ f(x) 



e{x - y) t 



p2t(y) 
\ f(v) 



<p(x)ip(y)dx dy 



p(x)p(y)e^i tk{ - xk+yk) e(x - y)ip(x)ip(y)dx dy. 



So, finally setting V(x) = \(V(x) + log /(a;)) yields 



T2n{t) = pf(m 2n ) 



(2n)\ Js 2n 



e Tr(-V(X) + j:?UX>) dx ^ 



The map O for the classical orthogonal polynomials at t = 0: Then, 
the matrix O, mapping orthonormal pt into skew-orthonormal polynomials 
qk, is given by a lower-triangular three-step relation: 



'c 2n 



V a 2n 



q2n+l(0,z) 



I02n 
C2n 



( c 2n 2n \ 

-C 2n -lP2n-l(0, Z) H -(52bi)p 2n {0, z) + C 2n P2n+l (0, z) 

\ fl 2n o / 

(6.2) 

where the a« and bi are the entries in the tridiagonal matrix defining the 
orthonormal polynomials, and the Ci 's are the entries of the skew-symmetric 
matrix M ■ 

In @, we showed that then Af is tridiagonal, at the same time as L, (see 
Appendix 2) 



b a 

a bi ai 

a\ b 2 



-Af 



c 
-c ci 

-ci 



(6.3) 



with the following precise entries: 



44 



A v M : Pfaff /matrix integrals March 8, 1999 #1 §6, p. 45 



Her mite : p(z) = e z * , a n ^ = b n = 0, c n = a n 



Laguerre : p(z) = e z z a I [0tQo) (z), a n _i = ^n(n + a), b n = 2n + a + 1, c r , 
a n /2 

Jacobi : p(z) = (1 - z) a {l + z) p I[- 1A] (z) 



An(n + a + (3)(n + a)(n + (3) 



1/2 



(2n + a + [3) 2 (2n + a + (3 + l)(2n + a + (3 - 1) / 
a 2 - (3 2 



Cr). (Ir 



(2n + a + (3)(2n + a + (3 - 2) 
'a + (3 



+ n + 1 



If the skew-symmetric matrix M has the tridiagonal form above, then one 
checks its inverse has the following form: 



/ 



■M- 1 = 




j_ 

CO 



gl 

C0C2 



C1C3 
C0C2C4 



C1C3C5 
C0C2C4C6 



_ J_ 





-ci 





-C1C3 





-C1C3C5 


co 


C0C2 


C0C2C4 


C0C2C4C6 





























_ J_ 





-C3 





-C3C5 




J_ 


C2 




C2C4 




C2C4C6 





















C2 






_ J_ 





















-C5 








J_ 


C4 




C4C6 





C3 
















C2C4 




C4 






_ J_ 
































C6 





C3C5 





C5 





J_ 







C2C4C6 




C4C6 




C6 





V 



(6-4) 

In order to find the matrix O, we must perform the skew-Borel decomposition 
of the matrix —LA 

-U = -Af- 1 = 0- 1 JO T -\ 

The recipe for doing so was given in theorem 4.1 (see also the important 
remark, following that theorem). It suffices to form the pfaffians (0.17), 
by appropriately bordering the matrix — A/" -1 , as in (0.17), with rows and 
columns of powers of z, yielding monic skew-orthogonal polynomials; we 
choose to call them r's, instead of the g's of theorem 4.1, with 0\{z) = 
r(z). They turn out to be the following simple polynomials, with l/f 2ri = 
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C C 2 C 4 • • • C 2 „-2 

1 C2nZ 2n 1 . o n 



r 2n{Z) = , = = -^=C 2n Z 



T 2nJ 2n+2 



c c 2 • • • c 2n wc 2 



n 



C 2n ^ — C 2n _i2; _ 1 

CqC 2 • • • C 2n ^/c^, 



r2n+l{Z) — i = = — — — \C-2nZ ~ C 2n -lZ )■ 

'T2 n T2n+2 



Then, also from appendix 1, in order to get O — > O in the correct form, we 
compute the skew-orthonormal polynomials f k , with Ox(z) = r(z): 



hn(z) = -Lr 2n (z) = J^z 2n 

V a 2n V fl 2n 



J2 h 

hn+iiz) = - — : r 2n (z) + y/a^r 2n+1 (z) 

V a 2n 



102 
C 2 



]n { J2n-l _|_ C 2n /V"^ j \ 2n , ^2n+l\ 

- -C 2n _iZ H {2_^ b i) Z + C 2n^ 

n \ «2n o / 

(6.5) 



From the coefficients of the polynomial f&, one reads off the transformation 
matrix from orthonormal to skew-orthonormal polynomials; it is given by the 
matrix O, such that Ox(z) = f(z). Therefore q(t, z) = 0(2t)p(2t, z) yields, 
after setting t = 0, 



<?2n(0,^) = J—P2n{0,z) 



a 2 



n 



q 2n +i{0,z) = J— 

V c 2n 



/ Q 2n \ 

-C 2n -lP2n-l(0, Z) H -(£ bi)P2n(0, z) + C 2n p 2n+1 (0, z) , 

V fl 2» / 

(6.6) 



confirming (6.2). 



7 Example 2: From Hermitean to symplectic 
matrix integrals 
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Proposition 7.1 The matrix transformation 

maps the Toda lattice r -functions with t- dependent weight 

Pt ( z )=e- y ^T^\ V' = g/f 

(Hermitean matrix integral) to the Pfaff lattice r -functions (Symplectic ma- 
trix integral), with t-dependent weight 

p t (z) := (p 2t (z)f(z))^ = e -^)-iog/(,)-2£~^) =; e -v(z)+EUz\ 
To be precise: 



p n (t, z) orthonormal polynomials in z for the inner-product 
(<f,^)t y = I <p(z)ip(z)p t (z)dz 



Toda lattice < 



fiijit) = (z\z J ) s t y , andm n = (/%)o<ij<n-i, 
r n (t) = detm n (t) = -( e W^^dI 



-N{2t) = 0-\2t)JO T -\2t) 
map 0(2t) such that ^ 0{2t) is lower-triangular 

0(2t)S(2t) G £ k 
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q n (t, z) — 0(2t)p n (2t, z) skew- orthonormal polynomials 

in z for the skew-inner-product (weight p t ), 



Pfaff lattice < 



sy 
2t 



p>ij(t) = {z\z 3 ) s t k , and m n = det(/^) <^<n-i 
f 2n (t) = pf(m 2n (t)) = I e^^^dX. 

Proof : Representing d/dx as an integral operator 



d 



dx ' 
compute 



x) = j^S(x-y)(p f (y)dy = - J -^-5(x-y)ip(y)dy = J^S'(x-y)(p(y)dy : 



= n t 



Pt dz 



'fp t , so that n t p(t,z) = Afp(t, z); 



remember Af from (5.2). Let it act on a function <f(x): 

Id 



/(*) x , 



X 



- y)\Jf(y)pt(y)v(y)dy. 



Then 



(<^>f = (<p,u*1;)% = (<p,n M il;)Z 



= I f P2t(x)<p(x)* -^p-5'(x-y)^ f(y)p 2t (y)i>(y)dxdy 



/ J R2 \] f(x)p 2t (x)ip(x)5'(x - y)^f f(y)p 2t (y)ip(y)dx dy 

R f (- 

/r 2 \ dx 



f(x)p2t(x)ip(x) J 5(x - y)\Jf(y)p2t(y)i>(y)dxdy 
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R 



' d_ 

dx 



\Jf(x)p 2t (x)(p(x) ^f(x)p 2t (x)^(x)dx 



+lJjf(x)p 2t (xMx) 



/ f(x)p 2t (x)ip(x) y-^y/ f(x)p2t(x)i/j{x)j dx 

= ~\ f {\/ f(x)f>2t(x)<p(x), ^ f{x)p 2t {x)i)(x)}dx 
= ~\ I W{x),^{x)}p 2 {x)e 2 ^ Uxl dx, 

using the notation in the statement of this proposition. Setting p(x) = e~ v( - x \ 
with V(x) = \{V{x) -log/) 

(x\x j ) sk = -If {x\x j }p 2 {x)e 2 ^ tixi dx 

2 JR. 

= -\ [ {x\x j }e- 2 ^-^ tiX ^dx, 

and so 

r 2n (t) = pf(m 2n (t)) = f e 2T ^v(*) + M dx . 

The map O 1 for the classical orthogonal polynomials at t — 0: Then, 
the matrix O, mapping orthonormal pu into skew-orthonormal polynomials 
qk, is given by a lower-triangular three-step relation : 



P2n(0, z) = -c 2n - 1 J^ L ^q 2n - 2 {Q, z) + ^a 2n c 2n q 2n (0, z) 

V C 2n -2 



p 2n +l(0,z) = -C 2 , 



uh^-q2n-2(°i z ) ~ (J2 b i)\ — 92n(0, z) + J—q 2n +i(0, z), 

(7.1) 



where the ai and hi are the entries in the tridiagonal matrix defining the or- 
thonormal polynomials, and the Ci the entries in the skew-symmetric matrix. 

In this case, we need to perform the following skew-Borel decomposition 
at t = 0, 

—U = -M = 0~ 1 JO T -\ 
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where J\f is the matrix (6.3). Here again, in order to find O, we use the 
recipe given in theorem 4.1, namely writing down the corresponding skew- 
orthogonal polynomials (0.17), but where the /j,^ are the entries of —U = 
—J\f: consider the pfaffians of the bordered matrices (0.17); they have leading 
term 



n-l 



T2n = II C 2i- 



Then one computes 
1 



r 2n - 

r 2n+l l 

with 
Setting 



p-=. — =51 

V T 2nl~ 2n+2 *=0 
1 



n /n—i—1 
2n-2i 



n-i 



II C 2j II C 2n-2j-l 



n-l 



' n—%—1 



n-i 



2n+l 



T~2n T 2n+2 



\J T 2nf 2n 



,2n-2i 



II C 2i II C 2n-2j-l 



i=l 



+2 = CoC 2 ...C2n-2V C 2n V T T 2 



CQ. 



D := diag(y r r 2 , y r r 2 , y r 2 r 4 , y r 2 r 4 , ...), 
the matrix O is the set of coefficients of the polynomials above, i.e., 



(7.2) 



O 



D 



-i 



( 1 


























1 




















Cl 





Co 

















c 2 








c 














Cl c 3 





c c 3 





CoC 2 











Ci c 4 





Co c 4 








Co c 2 








Ci c 3 c 5 





Co c 3 c 5 





Co c 2 c 5 





Co c 2 c 4 





Ci c 3 c 6 





Co c 3 c 6 





Co c 2 c 6 








Co c 2 c 4 



D^R 



(7.3) 



As before, in order to get the skew-symmetric polynomials in the right form, 
from the orthogonal ones, one needs to multiply to the left with the matrix 
E, defined in (8.2) in the appendix: 

6 = EO = ED^R, (7.4) 
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and so, 



0~ 



(7.5) 



it turns out the matrix O is complicated, but its inverse is simple. Namely, 
compute 



R- 1 



( 1 


_£l 

CO 
_£2 
CO 










1 













J_ 

CO 





C3 



CO C2 

C4 
CO C2 










J_ 

CO 












1 

CO C2 





CO C2 C4 

C6 
Co C2 C4 












1 

CO C2 















CO C2 C4 


















1 



\ 



Co C2 C4 



and 



V 



a 

-A) - 



o 



«2 
-fa ~ 



«4 

-Pa 





ad 



) 



(7.6) 



with a 2 n and /3 2ri as in (8.5). Carrying out the multiplication (7.5) leads to 
the matrix O" 1 , with a few non-zero bands, yielding the map (7.1). 
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8 Appendix 1: Free parameter in the skew- 
Borel decomposition 

If the Borel decomposition of —H = 0~ l JO T ~ l is given by a matrix O G Qk, 
with the diagonal part of O being 



/ 


O"0 










0"o 



(O)o 



\ 



0"2 








0"2 



a 4 








a 4 



■■ J 



then the new matrix 
/ 



O :-- 



1/ao 
(3 a 



l/a 2 
P2 a 2 



l/a 4 

j3 4 a 4 



V 



o 



O =: EO, 



J 



(8.2) 



with free parameters a 2n , fan, is a solution of the Borel decomposition — H 
0~ 1 JO T ~ 1 , as well. The diagonal part of O consists of 2 x 2 blocks 



l/«2n W a 2n 

$2n «2n/ V (T 2 ri 



0"2n/«2n 
/?2n0"2 n Q ; 2ri cr 2n 



Imposing the condition that 



q(z) = Op(z), with p k (z) = ^p^z 1 

i=0 

has the required form, i.e., the same leading term for q 2n and q 2n+ \ and no 



z 2n -term in q 2n +i, 



Q2n( z ) — <l2n,2nZ 2n + ' ' ' 

( ~\ — ~2n+l , ~2n-l , 

<?2n+H^J — q2n,2nZ + ?2n,2n-l^ + 



(8.3) 
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implies 



&2n 

&2; 



-P2n,2n ~ &2n&2nP2n+l,2n+l 



C2nf32nP2n,2n + 0'2n<^2nP2 



n+l,2n 







yielding, upon using the explicit form of the coefficients pu of the polynomials 
Pk, associated with three step relations (see next lemma), 



a 2n — 

fan 
«2n 



P2n,2r. 



P2n+l,2n+l 
P2n+l,2n 



P2n,2n 



0>2n 
Q>2n 



Hence 

So, if 
then 



Oi2n = \fa2^ and (3 2n = 
r(z) = Ox(z), 



1 2n 
&2n o 



(8.4) 



(8.5) 



/ 



r(z) := 6x{z) 



1/oq 



l/a 2 

fa «2 



l/a 4 

/?4 0:4 



and thus 



V 

f2n+l(z) 



r[z) 



Er(z), 



/ 0,2n 



-r 2 n(z) 



—=r 2n + Va2nr 2n+1 (z) 

\/(L2n 



(8.6) 
(8.7) 
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Lemma 8.1 A sequence of polynomials p n (z) = J27=oPmz l of degree n sat- 
isfying three-step recursion relation^ 

Zp n = O-n-lPn-l + KPn + O-nPn+l, U = 0, 1, . . . , (8.8) 

has the form 

p n+1 (z)=^(z n+1 -(±b i )z« + ..). 
Proof : Equating the z n+l and z n coefficients of (8.8) divided by p n>n yields 

Pn+l,n+l _ I 
Pn,n <^n 

and 

Pn,n—1 Pn+l,n , 7 

— ■ = a n h o n - 

Pn,n Pn,n 

Combining both equations leads to 

Pn+l,n Pn,n—1 _ , 

Pn,n Pn—l,n—l 

yielding 



a 



Pn+l,n 
Pn,n 



n 



n 



^2bi, using a_i = 0. 



o 



9 Appendix 2: Simultaneous (skew) - sym- 
metrization of L and J\f . 

For the classical polynomials, the matrices L and M can be simultaneously 
symmetrized and skew- symmetrized. 

We sketch the proof of this statement, which has been established by us 
in H. Given the monic orthogonal polynomials p n with respect to the weight 
p, with p'/p = —g/f, we have that the operators z and 



[fd FT f d ,f'-9 
n = \ — rvfP = f-r + 



p dz^ dz 2 



2 with a_i = 0. 
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acting on the polynomials p n 's have the following form: 

Zp n = d^Pn-i + b n p n + p n+1 
Ilp n = ... -JnPn+l, 



(9.1) 



in view of the fact that for the classical orthogonal polynomials^, 
Hermite: n = 4- — z 

dz 

Laguerre: n = z-^ — \{z — a — 1) 
Jacobi: n = (1 - z 2 )£ - ±((a + (3 + 2)z + (a - (3)). 
For the orthonormal polynomials, the matrices L and —A/" are symmetric 
and skew-symmetric respectively. Therefore the right hand side of these 
expressions must have the form: 



a„_lPn-l + KPn+Pn+1 
2 



\\p n = 0-_ 1 7n-lPn-l -InPn+l- 

Therefore, upon rescaling the p n 's, to make them orthonormal, we have 
zp n = (Lp) n = a n _ x p n _i + b n p n + a n p n+1 

lip n = (J\fp) n = a n „i7 n _ip n _i - OnlnPn+l, 

from which it follows that 

c 
-c ci 

-ex 



-Af 



with c„ = -a n j n 



where — 7„ is the leading term in the expression (9.1). 



10 Appendix 3: Proof of Lemma 3.4 

For future use, consider the first order differential operators 



00 z~ j d 

v(t,z) = Y,—— and B ( z ) 



dz 



having the property 

B(z)e-^f(t) = B(z)f(t - [z' 1 ]) = 0. 



10.1) 



10.2) 



2.1 



with respective weights p — e 2 , p — e z z a , p = (1 — z) a (l + z) 13 . 
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Lemma 10.1 Consider an arbitrary function (p(t, z) depending on t £ C°° ; 
z G C, having the asymptotics <f(t,z) = 1 + 0(~) for z f oo and satisfying 
the functional relation 

^-^Vi) = ^-[^r 1 U) > teC -, zee. (10.3) 

(p{t,Zi) ¥{t,z 2 ) 
Then there exists a function r(t) such that 

¥>M = • (10-4) 

Proof : Applying Si := B(z\) to the logarithm of (10.3) and using (10.1) 
and (10.3) yields 

(e-^ -l)B 1 \ogi P (t;z 1 ) = -B x log<p(t, z 2 ) 

oo 



j= i ot j 

which, upon setting 

fj(t) = Res zl = 00 z{B 1 log(p(t,z 2 ), 

yields termwise in zi, 

(e-iiCa) _ l)f 3 (t) = -JLl 0g(p (t, Z2 ). (10.5) 

Acting with on the latter expression and with ^- on the same expression 



with j replaced by i, and subtracting^, one finds 

'dfi dfA 



(e-' n{z2) - 1) 

yielding 



v dtj dU 



df % df 



j 



dtj dti 

the constant vanishes, because fr 1 never contains constant terms. 



24 It is obvious that 



0. 
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Therefore there exists a function \ogr(t 1 ,t 2 , ...) such that 
d 

and hence, using (10.5) 

Alog^ )Z ).(e-^-l)Alogr 
or, what is the same, 

J-(log V 9-(e- r '-l)logr)=0, 



from which it follows that 



00 bi 



\og V -(e^-l)\ogr = -J2-z- t 

l 1 



is, at worst, a holomorphic series in z 1 with constant coefficients, which we 
call —bill. Hence 

<p{t,z) = 1 



T(t) 

rit-iz- 1 ^^-^) 



i.e. 



ip(t,z 



f(t-[z-i}) 



where 



f(t) 
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